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Abstract 

The influence of strong and electroweak penguin amplitudes m B / B tt+tt" is investigated in 
connection with the determination of the unitarity triangle angle a of the CKM matrix. A rela- 
tion between the observable asymmetry, the angle a, and the penguin amplitude is established. 
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A model calculation of the penguin amplitude shows that the CP asymmetry m. B^ ^ -k 
decays is only mildly influenced by the penguin amplitudes. Experimental limits on pure pen- 
guin and penguin dominated processes are consistent with the model. This information also 
suggests in a rather model independent way that penguin amphtudes will not be a serious 



5^ \ complicating factor in the determination of a from the 7r"*"7r time dependent asymmetry. 
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1 Introduction 



It is expected that B decays will show large CP- violating effects, characterized by non- vanishing 
values of the angles a, (3 and 7 in the unitarity triangle [0]. One of the best ways to detect this 
CP violation is to measure an asymmetry between and decays into a CP eigenstate. If 
only one weak amplitude contributes to the decay, the phase in the elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix can be extracted without uncertainties due to unknown 
hadronic matrix elements. Thus sin 2a, sin 2/5 and sin 27 can in principle be measured in 
5°, B^ TT+TT", J/ipKs and Bs, Bg — > p'^Ks decays, respectively. Unfortunately the situation 
is more complicated. In all of the above cases, in addition to the tree contribution there are 
amplitudes due to strong and electroweak penguin diagrams. In the case of the J/ipKs final 
state the weak phase of the penguin term is the same as that of the tree contribution. Thus 
there is no uncertainty for determining sin 2/5 from the CP asymmetry. 

For 5°, B^ 71^71^ the weak phases of the tree and penguin contributions are different 
causing hadronic uncertainties in the interpretation of an otherwise clean experiment. How- 
ever, by measuring also the rates of B^ tt^tt^, B^ tt^tt^ and their charge conjugate decays 
one can isolate the amplitudes contributing to final states with isospin and 2 and thereby 
determine a . This construction, however, relies on the fact that electoweak penguin contri- 
butions do not exist, since they contribute to both isospins and not only to / = as the strong 
penguin terms 0. Although these weak penguin terms are expected to be small compared to 
the tree amplitudes 0, so that the Gronau-London construction should be possible, there is 
still the problem that the partial rates of the decays B^, B^ — > vr^vr^ are at least an order of 
magnitude smaller than for the other 27r final states 0. In addition, because of two neutral 
pions in the final state, these decays are very difficult to measure accurately. So if this program 
can not be carried out the error of sin 2a is of the order of \P/T\, where P (T) represents the 
penguin (tree) contribution to B^ tt+tt^. In this connection DeJongh and Sphicas studied 
the behavior of the asymmetry based on a general parameterization of the penguin magnitude 
and phase 

Recently, two of us p calculated the effect of strong and electroweak penguins in all B"^'^ — *• 
TTTT, ttK and KK decays using specific dynamical models for the hadronic matrix elements. 
Concerning the asymmetry between r{B^ tt+tt") and r{B^ tt+tt^) (Acp) it turned 
out that the effect of electroweak penguins was indeed small, of the order of 2%, and that 
strong penguin amplitudes changed the asymmetry by less than 20% as compared to the tree 
value. These results were fairly independent of the specific models employed for calculating the 
hadronic matrix elements. Since the parameters in the time dependent asymmetry are obtained 
from ratios of the weak transition matrix elements, it is clear that they are much less model 
dependent than, for example, the branching ratio. Of course, this rather moderate change of 
Acp for B'^ — > TT+TT" depends on P/T which was determined by the model calculations. Due 
to the way the results in were presented only one particular set of CKM parameter values, 
namely p = —0.12, rj = 0.34 was assumed. Although this is the preferred value obtained in 
the analysis of 0] in their so called "combined fit" it is certainly not the only possible set 
following from their analysis. From CP violation in the — system it is known that 
rj ^ 0. Nevertheless, for both p and rj, only very loose bounds exist which translate into similar 
loose bounds on the triangle phases a, (3 and 7 

10° < a < 150°, 5° <P < 45°, 20° < 7 < 165°. (1) 

From some predictive SUSY GUT models on fermion masses and mixings, a was found to be 
large [|, ^. For example, the model proposed in predicted: a = 86°, /? = 22°, 7 = 72°. 
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It is clear that the change in Acp for B^/B^ tt+tt^ due to penguin contributions depends 
not only on P/T but depends also on the particular set chosen for p and t]. However, assuming 
fixed p and t] values in |^ was an unnecessary limitation. Of course, it would be easy to repeat 
the calculation of for any other set of p, r] inside the bounds of (1). This would give us a large 
array of numbers for Acp- Instead we follow in this note a different route, which is particularly 
simple when we neglect the electroweak penguin terms and use some approximation for the 
strong penguins. We express the main contribution to the asymmetry parameter, fle+e'? which 
is the coefficient of the sin(Amt) term in Acp (see below) in terms of the tree and penguin 
amplitudes and their relative phase. Then a^+e' depends only on a. This gives us a clear insight 
into the dependence on \P/T\ and on the strong phase and allows us to derive upper limits on 
the change of cte+e' by including information from other decay channels which depend on the 
penguin contributions more strongly than the decay into vr+vr". 

The outline of the other sections is as follows. In section 2 we give the formulas of the 
asymmetry, from which we start and derive the formula for the change due to the penguin 
terms. In section 3 we present our results and discuss their relevance. 



2 CP-violating Observables in tt+tt 

In this section we establish a relation between the CP- violating observables in B^ /B^ tt+tt", 
the angle a of the CKM matrix, and an auxiliary variable ao involving the ratio of penguin to 
tree amplitude and the strong interaction phase difference between the tree and the penguin 
amplitude. Applying the general analysis on rephase-invariant CP-violating observables given 
in ref. |]T^ for the B-system, and expressing the two physical mass eigenstates Bl and Bh as 

Bl = p\B° > +q\B^ >, Bh = p\B° > -q\B^ > (2) 

with the decay amplitudes of B^ tt+tt^ and 5° tt+tt^ written as 

g =< 7i+7i-\Heff\B^ >= ATe'^^+'^^ + Apd'^''^'^^ = h, (3) 
h =< 7r+7r-|/7e//|i?° >= Atc'''^-^^'^-^ + Ape'"'^''^'^'' = g (4) 

the time-evolution of states with initially pure B° and B^ are found to be 



T{B%t) TT+TT-) oc T^^^^^e-r*[(l + a,a,0 cosh(Art) 
+ (1 + Qee') sinh(Art) + (oe + a^i) cos(Amt) + a^+e' sin(Amt)] (5) 

T{B%t) TT+TT-) oc ^i^^^e-r*[(l + a,a,0 cosh(Art) 
-|- (1 + a^e') sinh(Art) — (a^ + a^/) cos(Amt) — a^+e' sin(Amt) ] (6) 

where a^, a^', a^+e' and a^^' are rephase-invariant observables and defined as follows 

1 - _ 2ReeB _ 1 - \h/g\'^ _ 2Ree's 



~AIm{qh/pg) 2/meB(l — < 


1- 

~b\ 


+ 2/me^ 


(1- 






(l + \q/p\2)(l + \h/g\^) (1 + 




m + 







(7) 



ARejqh/pg) _ AlmeB Ime'^ - 2{\eB\^ + \e'^\^) 

{l + \q/p\2)(l + \h/g\^) (l + |e^|2)(i + |,/^|2) 
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with = (1 — q/p)/ (1 + q/p) and e'^ = {1 — h/ g)l{l + h/g). In the B system, since <C 1, 
|Ar| ^ I Am I and |Ar/r| ^ 1 , the time-dependent asymmetry Acp(t) can be simply written 

Y(B° f) -T(B° ^ f) 
Acp{t) = rj^o^/) + r(Bo^/) - ^°'(^"'^) + «^+^' sin(Amt) (8) 



The CP-violating phase is related to the observables via 



sin(2(0M + 0^)) = , ^ ^ (9) 



where the phase (pM and are defined by 

P V 9 9 

The tree amplitude is proportional to f „ whereas the penguin amplitude depends in general 
on Vu and fc, where Vu = VubV*^, Vc = VcbV*^ and Vt = VtbV^^. It is well known that when the 
difference of the u and c contributions in the qq intermediate states can be neglected, the penguin 
amplitude can be expressed in terms of Vt alone. A general analysis of these considerations has 
been carried through by Buras and Fleischer |jTl]]. In the model to be considered in the next 
section this is only violated by the additional 0{as) and 0{q) corrections in the short distance 
coefficients @. In this approximation, for — > tt+tt" decay, we have 

h = Vuf + vtP 

= \vu\Te-'^ + \vt\Pe'^ (11) 

Then we have 



M 


= p 




= 7 


h 


= -/5 



(12) 

where a, (3 and 7 are three angles of the unitarity triangle of the CKM matrix and a = 'K — f]—^. 
Factoring the phase (pT of the tree contribution out we introduce the phase shift due to the 
penguins, ao, defined by: 

0A = 0T - ao (13) 



As a result we can write 



sin(2(a + ao)) ^ a^+e' (14) 



where ao is given by 

tan 2ao 



2(4f)sinA0 cosA,5+(4g)^sin(2A0) 
l + 2(^)cosA0 cosA<5 + (^)2cos(2A</.) 

2{^)sma cos5-{^)'sm{2a) 
l-2(^)cosa cos5 + (^)2cos(2a) 



(15) 
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with A0 = 0r — 0p = 7r — a the weak phase difference and AS = 5t — Sp = S the strong phase 
difference between tree and penguin diagrams. 

When the strong phase S is zero, this equation simphfies to: 

Ap sin a 

tanao = —- (16) 

(14) is a relation between the asymmetry a^+e', the unitarity triangle angle a, and the 
penguin complication represented by ao- (15) shows how the angle ao depends on the size of 
the penguin, the strong phase 6 and the unitarity angle a itself. To determine a from the 
B^/B^ TT~^7T~ time dependent asymmetry, must be calculated from a model or estimated 
from some other process. This is the subject of the next section. 



3 Extraction of a 

The three angles a, (3 and 7 are defined as 

"="'-n"K^J' ''""'■n-vyaj' ''="n"^J' '''' 

They are related to the Wolfenstein parameters p and r] as follows: 

tana = r, tan/3 = , tan7 = - (18) 

r-p(i-p) 1-P P 

From the present experimental data on |Kjfe/Kb|; one has [0 



p2 + ^2 = 0.36±0.08 = (19) 

It is clear from (14) that to extract a from experiment, one has to know ao. As shown in (15), 
ao depends on a (or rj and p) as well as the ratio of the penguin amplitude to the tree 
amplitude At- We can separate the CKM matrix elements and pure hadronic matrix elements 
in the amplitudes Ap and A^ in the following way. 

At = T, Ap = \vt\P (20) 

(20) is correct when we neglect the difference between the u and c quark contributions to the 
penguin amplitude, discussed further below. In this same approximation, the strong phase S is 
zero in our model; then in this limit (15) can be further simplified to: 

'^''''' = B--{B^-p)PlT Kf) ^^'^ 

In general, calculating the strong phase 6 is difficult due to unknown nonperturbative effects. 
In 1^ the strong phases derive from absorptive parts of the qq intermediate states in the strong 
penguin contributions which are estimated perturbatively using recently developed next to 
leading log formalism following the method pioneered by [0. In it is found that 



Sp ~ 9.5°, 6 ~ -9.5°, 6t = (22) 
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For a general consideration, we take 5 as a free parameter. The ratio (P/T) is purely determined 
by the hadronic matrix elements. In the operator product expansion approach, the hadronic 
matrix elements are products of short-distance parts, i.e. Wilson coefficients, evaluated by 
perturbative QCD, and a badly known long-distance part. In a factorization approximation 
the long distance hadronic matrix elements are themselves products of current form factors and 
decay (coupling) constants. For the tt+tt^ decay, it is easy to see that the ratio {P/T) 

is almost independent of the uncertainties in the long-distance modeling because differences in 
the approaches cancel in the ratio. Therefore, can be well determined in a rather model 
independent way and given by coefficients a, which have been defined in in terms of the 
effective Wilson coefficients c/^ : 

= 4'' + ^cf ^ (23) 

where (i,j) is any of the pairs (1,2), (3,4), (4,5), (6,7), (7,8) and (9,10) and N is the number of 
colors. (See |^ for further details.) The second term in (23) arises from the Fierz rearrangements 
in connection with the factorization contributions. In two models were considered, N = oo 
and = 2, to account for possible non-factorizable contributions. Then as one can see from 
Tab. la,b of the ratio \P/T\ becomes: 

P |a4 + Qio + (q6 + a8)-R[7r~,7r+]| 

As we see from (24) this ratio does not depend on the current form factors and decay constants. 
It depends only on the factorization hypothesis and on the effective short distance coefficients. 
Furthermore it is found that the ratio is not sensitive to the effective color number N in (23). 
The reason for the simple structure of \P/T\ as given by (24) is that for 7r"'"7r~ states there is 
only one way to factorize the transition matrix element. 

The tree and penguin amplitude for tt+tt^ were evaluated in [^] for various hypotheses 

concerning the 0{as) corrections in the c^-^'^ (which include the absorptive parts) and the 
influence of the electroweak penguins. From these results we can calculate the amplitudes T 
and P with the CKM phases factored out. Since the current matrix elements cancel in the ratio 
this is equivalent to evaluating P/T in terms of the coefficients Oj. The results are displayed in 
Tab. 1 for N=2 and N=cx3. The notation Py refers to penguin amplitudes arising from x = st 
(strong) or X = ew (electroweak) penguins, and y = u or y = c parts of the weak Hamiltonian 
with and without 0{as) terms in c^-^'^ , where the absorptive parts are contained in the 0{as) 
corrections. In the following we shall use these results in order to calculate for various 
assumptions concerning 0{as) terms in the c^^^ or electroweak penguin effects. The simplest 
case is to use the tree and strong penguin amplitudes with 0{as) corrections neglected, i.e. the 
penguin amplitudes of column 3 and 4 in Tab. 1 which results in 5 = 0. The relation between 
these amplitudes and the penguin introduced above is P = — |(Pu* + Pf )• For this case ^ in 
(24) is 0.05 for both N=2 and N= oo. For fixed B (21) determines Oq as a function of a using 
(18). 

In Fig. 1 ao is plotted as a function of a with B = (p^ +77^)2 = 0.28, 0.36, 0.44, respectively, 
for 5 = 0. From Fig. 1 it is apparent that the penguin shift decreases with increasing B. The 
maximum of ao as a function of a occurs near p = 0. In the plot p varies with a starting from 
p < to p > with increasing a. Fig. 1 contains our main result. As one can see the penguin 
shift, ao, is small compared to a. ^ is largest at small a and decreases monotonically with 
increasing a to zero. 

To see the influence of a strong phase 6 we have repeated the calculation of a now using (15) 



6 



where we still use the same ratio of P/T from column 3 and 4 in Tab. 1. In Fig. 2. the angle 
is plotted as a function of a for B=0.36 and four values of the strong phase, 5 = 0°, 10°, 40°, 90°. 
The maximal distortion occurs for 5 = 0° as is evident already from (15). The behavior for 
90° < S < 180° is approximately a reflection ao(5) = —ao{7i — 6). ao{6) is an even function of 
5. 

The relation between a and the measured angle is a = an — cto: which can be read off 
from Fig. 3. This plot is for 6 = 10° and B= 0.28, 0.36, 0.44. As we can see, a as a function 
of ajvf = a + aQ which comes from the measurement of a^+e' is rather independent of B and 
differs only slightly from the straight line of slope unity it would be, if there were no penguin 
contributions. 

The results for presented so far are for Wilson coefficients in which 0{as) and electroweak 
contributions have been neglected. The 0{as) terms yield deviations from (20) and additional 
absorptive contributions which generate the phase in (22). In Tab. 1 the penguin terms with 
the 0{as) terms included are given in the fifth and sixth columns as P^* and P^^. As we can 
see the 0{as) terms change the penguin terms up to 24% and produce the deviation 

pst pst 

= 0.13 

pst I pst 
u ' c 

in the real parts. In addition there is an imaginary part of the same order as the real part 
(these numbers are for the N=2 case). Of course the resulting shift in depends on a or 
equivalently on the value of p. Instead of calculating ao as a function of a (or p) we quote only 
results where a nearly had its maximum, i.e., p = 0, at a = 70°. By calculating a[ and ag+g/ 
directly we can use (14) to extract ao for this value of a. The major effect is due to the strong 
penguin amplitude itself, without as or electroweak corrections, which shifts a by ao = 8.0°. 
The Us and electroweak corrections shift a by an additional Aa = 1.9° and 0.6° respectively, 
to a total shift of ao = 10.5°. 

The values of the strong phase may be extracted from the direct CP-violating parameter 
Oe'. Hence in principle one can determine 5 without recourse to a model. On general grounds, 
however, one would expect 5 to be small. 

It is clear that the shift Uq due to the penguin effects is small. For all a the relative shift 
ao/a is less than 30% and decreases from this value with increasing a. This result depends on 
our model of the penguin amplitudes, in particular on P/T. In principle one need not rely on 
the model but rather obtain \P\ from pure penguin or penguin dominated processes, i.e., K^n~ 
(pure penguin) and K^n^ (penguin dominated). Unfortunately there are only experimental 
upper bounds on the branching ratios of these decays. In reference we have calculated the 
branching ratios of these decays for the special choice p = —0.12, r] = 0.34. Since K^7t~ depends 
only on \Vts\ which is well known, we can obtain upper limits on \P\. Taking for example the 
N=2 case, when we compare with the experimental limit |13|, BR{B- A'^tt") < 4.8 x 10-^ 
we find that \PexplP\ < 2.2. Keeping T fixed, such a large value of the penguin amplitude given 
by the upper limit of \Pexp\ would increase from 8° to 17° for the maximal shift at p = 0. Here 
we make the assumption that the penguin amplitudes of Ktt and vrvr final states are related, 
i.e. a larger P in Kii would mean a larger P in vrvr. This can be justified with SU(3) symmetry 
arguments. SU(3) symmetry breaking effects are indeed moderate in our model calculations. 

The penguin dominated decay channel — > K^Ti~ gives us a much better limit as ad- 
vocated by Silva and Wolfenstein [|14|. The experimental limit on this branching ratio is 



1.7 X 10 ^; from our previous work we know that this decay is dominated by the penguin 
amplitude (since the tree amphtude is Cabibbo suppressed) in the ratio 4:1 in the ampli- 
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tude. This gives about \Pexp/P\ < 1-4 leading to an even smaller shift of to ~ 11° com- 
pared to the 8° calculated above. Other experimental limits on the branching ratios relevant 
for comparing with our model are BR{ti~tt^) < 1.6 x 10~^, BR{'k^t[~) < 2.0 x 10~^ and 
BR{'K^TT~ + K^TT~) = (1.8 ±0.6) X 10~^ |]13|. Our model obeys these constraints; in particular 
we obtained BR{t[^tt~ + K'^tx^) = 2.15 x 10^^ for p = 0,1] = B which agrees perfectly with 
the measured value. This shows that the tree and penguins can not be too far from the values 
in our model. 

In conclusion, we find that the penguin distortion in the determination of a from the tt+tt" 
asymmetry is not a real obstacle provided a is not too small. Even if the penguin amplitude is 
taken from the upper limit of the pure penguin dominated process the penguin distortion on a 
remains below ~ 25% at large a. Improved experimental work on exclusive charmless hadronic 
B decays will even more sharply constrain the size of penguin amplitudes and in turn limit the 
shift aQ. 
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Table Caption 

Tab. 1: Reduced amplitudes for — > tt+tt^. CKM phases have been removed. These ampli- 
tudes are based on a next to leading log corrected weak Hamiltonian with as and a corrections, 
assuming factorization of the hadronic matrix elements. Numbers in the parenthesis are the 
real and imaginary part of the amplitude. For further details, see reference 0. 
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Tab.l 





Tree 


Reduced Amplitudes B'^ vr+vr 
Strong and EW Penguins T, P**, P^"' with () or without (") corrected 

NLL QCD Coefficients 


N 


T 


pst" 


pst" 
c 


pst 
u 


pst 


pew 


pew 

c 


oo 
2 


2.84 
2.44 


(-0.141,0) 
(-0.115,0) 


(-0.141,0) 
(-0.115,0) 


(-0.134,-0.0548) 
(-0.109,-0.0457) 


(-0.175,-0.0291) 
(-0.143,-0.0243) 


(0.0052,0) 
(0.0052,0) 


(0.0052.0) 
(0.0052,0) 
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Figure Caption 

Fig. 1. tto as a function of a with 5 = 0°. The curves correspond to B = 0.28 (dotted), 
0.36 (sohd), and 0.44 (dashed). All angles are in degrees. 

Fig. 2. aQ as a function of a with B = 0.36. The curves correspond to the strong phase 
(5 = 0° (dashed), 10° (solid), 40° (dotted), and 90° (dot-dashed). All angles are in degrees. 

Fig. 3. a = aM — ao as a function of the measured q;m with 6 = 10°. The curves corre- 
spond to B — 0.28 (dotted), 0.36 (solid), and 0.44 (dashed). All angles are in degrees. 
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